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Abstract 

We consider the low-temperature behavior of NMR relaxation rates T\ and 
T2G in the spin-1/2 antiferromagnetic chain. We find that T\ oc const and 
T2G oc y/T, with logarithmic corrections. We determine both constants and 
logarithmic terms by matching perturbative RG with exact results, so that the 
final expressions for the relaxation rates do not contain any free parameters. 
Our theoretical results are in excellent agreement with NMR experiments in 
Sr2Cu03. 
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There has been significant interest in the low- dimensional quantum spin systems and 
quantum critical phenomena in recent years due to the discovery of high-temperature super- 
conductivity in layered cuprates. Nuclear magnetic resonance (NMR) has been a powerful 
tool in studying the spin dynamics of the cuprate compounds jl]. The measurements of 
the longitudinal relaxation rate 1/Ti, the spin echo decay rate 1/T 2 g, and the Knight shift 
over a wide temperature range define the spectrum of antiferromagnetic fluctuations in the 
high-T c cuprates || . These NMR experiments have been the basis of the phenomenological 
spin-fluctuation model ||, which exhibits quantum critical behavior [p]], because the system 
always remains close to the antiferromagnetic instability. 

Quantum spin- 1/2 spin chain is another example of a quantum critical system with novel 
properties. The physical properties of this system are well understood theoretically. The 
XXZ spin-1/2 quantum spin chain is described by the Hamiltonian: 

H = JJ2[S*S? +1 + S!S? +1 + yS*S* +1 ). (1) 

i 

This model exhibits quantum critical behavior for —1 < 7 < 1, with asymptotic correlation 
functions vanishing with distance as a power law. Inverse temperature 1/T acts essentially as 
a finite system size, so that the power law for correllation function crosses over to exponential 
decay at distances r > £t — hc/T For I7I > 1 the spin-1/2 chain is equivalent to a model 
of massive free fermions, and the correlation functions decay exponentially with distance. 
Exactly at the Heisenberg point, 7 = 1, logarithmic corrections to the correlation functions 
appear as a result of the presence of the leading marginally irrelevant operator. 

Logarithmic corrections in the SU(2)-invariant models are well known. For example, a 
calculation of the logarithmic corrections up to two loops was done for the fermion model 
with backward scattering [[|, the sine Gordon model ||, and the SU(2) Gross-Neveu model 
||. Logarithms appear in every physical property of the spin-1/2 Heisenberg chain. For 
example, 1/ ln(if ) dependence of magnetization on magnetic field was found in Ref. [[!]] and 
in Refs. The temperature-dependent l/ln(T) corrections to the bulk spin suscepti- 

bility 0. Finite -size scaling corrections to the staggered spin susceptibility were found in 



Ref. [p| . In what follows we determine the logarithmic behavior of the nuclear magnetic 
resonance (NMR) relaxation rates 7\ and T 2 q. 

Nuclear spins probe local spin environment. The Knight shift provides a measure of the 
uniform magnetic susceptibility at a particular nuclear site, while the experiments on the 
spin-lattice relaxation rate yield information on the imaginary and real parts of the dynamic 
spin susceptibility x(q, u;). The analysis of the NMR experiments begins with the magnetic 
hyperfme hamiltonian, which couples nuclear spins and conduction electron spins: 

H HF = ^haSja, (2) 

I is the nuclear spin, S is the electron spin, a enumerates spin projections for sites i and j 
. The following expressions can be obtained for Ti and T 2 g'- 



j_ = 2^T rdq Im X (q,u ) 

T x h 2 J 2tt ±w u V ; 
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Here ^4||(g) and A±(q) are the hyperfme couplings parallel and perpendicular to the easy axis 
of the crystal, c^o is the nuclear resonance frequency, which is much smaller than any other 
electron energy scale. The magnetic field is directed along the c-axis. The q-dependence is 
smooth and arises from appropriate form factors. The susceptibility \ should, in principle, 
include contributions from both the uniform and staggered spin fluctuations. However, 



simple power counting |13| shows that the staggered component is dominant at small T. 
Indeed, the contribution of the uniform component scales as 1/Tj oc T, 1/T 2 g oc T°, while 
for the staggered component 1/Xi oc T°, 1/T 2 g oc T~ x l 2 . 

For the purpose of comparison of theory with experiment, it is convenient to define 
normalized dimensionless NMR relaxation rates: 

hJ (5) 
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The calculation of the NMR relaxation rates are based on the continuum limit bozonised 
approximation to the Heisenberg model. The Hamiltonian density can be written as: 

H = H - 2ngjJ R . (7) 

Here Hq is the Hamiltonian density for a free boson of compactification radius R = 1/ V2n, 
Jl and Jr are the left and right moving currents. The coupling constant obeys the renor- 
malization group equations known for the Kosterlitz-Thouless or the Kondo problem with 
ferromagnetic interactions [fbl: 



P = dg/dlnA=-g 2 -±g 3 , (8) 



as the energy cutoff A varies. 

The low-temperature asymptotics for NMR relaxation rates are determined by the stag- 



gered component of the dynamic spin correlator, which has the following general form |15 



Xs (r,t,T,g ) = DZ(g(T))F[g(T),rT/c,tT}. (9) 

Here Z(g(T)) is a cutoff-dependent renormalization factor, C is a nonuniversal constant. 
This scaling form can be easily obtained from the Callan-Symanzik Renormalization Group 
equations: 

[-d/d In T + P(g)d/dg + 2rf{g)] X ,(r, rT/c, tT, g) = 0, (10) 
where (3{g) is the beta function for the coupling constant g in Eq.(^): 

-f3(g)=g 2 -g 3 /2 + 0(g±) (11) 



dlnT 

and 7(g) ~ 1/2 — g/4 + 0(g 2 ) is the anomalous dimension. In Eq. ( |T0| ) the T-derivative 
acts only on the first argument of Xs] r T and tT are held fixed. The solution of Eq.(|T0"D can 
then be written in the form Eq. (|9]) with 

D K exp ( 2 f m dg ) (12) 



where go = g(A) is the "bare" coupling - the coupling at the energy cutoff scale A. 

The leading behavior of Z(g) is easily seen from perturbative expressions for (3(g) and 
7(0) given above fllB] , Z(g) oc \j yfg. Integrating Eq.(|9|) gives 



- oc Z(g(T))(l + a l9 (T) + a 2 g(T) 2 + • • •) 
1 1 

-L oc ^221(1 + M(T) + MT) 2 + -..). 

-t2G V -I 



(14) 



The non-universal constants dj, fej, the cutoff A and common factor D in these expressions 
can be fixed using exact Bethe Ansatz results on the correlation functions. The value of 
D = 1/(2tt) 3 / 2 was determined in Ref. 0>IH • 



A complete calculation of the NMR relaxation rates gives: 
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Here D = l/(2vr) 3 / 2 is the non-universal amplitude, C ~ 0.5772157 is the Euler's constant, 
while the integrals Jo and I\ are given by 



71.2766 




-259.94, 



(17) 



A J V 4 

where ty(x) is digamma function. The 1/ ln(A/T) term is incorporated in the cutoff as in Ref 
P~5|1 . Thus, up to terms 0(1/ ln 2 (A/T)) the temperature dependence for 1/Tj or v / T/T 2 g is 



actually given by the square root of the log and loglog terms in the numerator of Eq.(|l7]). 
The cutoff parameters A can be determined using exact methods: 

A = 2 v / 2^e c+1 J ~ 24.27J (18) 

Ai= V^ e - h/ 2I OJ ^ 5 27 j (1Q) 



The ratio of the relaxation rates, however, is only weakly temperature dependent, 

( T *>\ - LSSofl + S) (20) 



\TiVTj ' V ln(A/T)/ 

\ 1 v / norm \ \ I / / 

In summary, I have presented theoretical low temperature results for NMR rates T\ and 
T2G-The temperature dependence of 1/T% is logarithmic, which is expected in the presence 
of the leading irrelevant operator. The main improvement over past theories is that the 
amplitude and the cutoffs were also calculated, so our expressions do not contain any free 
parameters. This has been a matter of confusion in the past [ 19fl . In particular, cutoffs 
have been fitted and used together with 1 /log terms which redefine them. Presumably, 
this explains the difference of our cutoff A with fitted value of 4. 5 J used in Ref. [[HJ. Our 
theoretical results agree reasonably well with experimental data of Takigawa et al. j[9] on 
Sr 2 Cu0 3 , as shown in Fig.l and Fig. 2. 

This work was supported by the National High Magnetic Field Laboratory through NSF 
cooperative agreement No. DMR-9527035 and the State of Florida. 
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FIG. 1. NMR 1/Ti for magnetic fields along different crystal axes in S^CuOa from Takigawa 
et al. [ 19 1 compared to our theoretical expression with no adjustable parameters. 
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FIG. 2. NMR T l / 2 /T 2G from Taki gawa et al. |l9| vs our low-temperature result. The data 
shows different assumptions for the value of J (fitted from the spin susceptibility data). The 
crosses are the data for J = 2850-fT not corrected for the I z fluctuations. 
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